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We present a method for calculating the radiative tranfer on a protoplanetary 
^ ■ disk perturbed by a protoplanet. We apply this method to determine the effect on 

the temperature structure within the photosphere of a passive circumstellar disk 
in the vicinity of a small protoplanet of up to 20 Earth masses. The gravitational 
o potential of a protoplanet induces a compression of the disk material near it, 

q ■ resulting in a decrement in the density at the disk's surface. Thus, an isodensity 

+3 . contour at the height of the photosphere takes on the shape of a well. When 

such a well is illuminated by stellar irradiation at grazing incidence, it results 
in cooling in a shadowed region and heating in an exposed region. For typical 
^ ■ stellar and disk parameters relevant to the epoch of planet formation, we find 

that the temperature variation due to a protoplanet at 1 AU separation from its 
parent star is about 4% (5K) for a planet of 1 Earth mass, about 14% (19 K) for 
planet of 10 Earth masses, and about 18% (25 K) for planet of 20 Earth masses, 
We conclude that even such relatively small protoplanets can induce temperature 
variations in a passive disk. Therefore, many of the processes involved in planet 
formation should not be modeled with a locally isothermal equation of state. 
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1. Introduction 

Planetary systems are formed from rotating protoplanetary disks, which are often mod- 
eled in steady-state, in vertical hydrostatic equilibrium, with gas and dust fully mixed and 
thermally coupled (Kenyon & Hartmann 1987). Such disks describe very well the observed 
properties of T Tauri disks of age ~1 Myr and typical mass accretion rates of < lO _8 M yr _1 
(Hartmann et al. 1998; D'Alessio et al. 1998). They are passive disks in the sense modeled by 
e.g. Chiang & Goldreich (1997), where the main source of photospheric heating is irradiation 
from the parent star, although viscous heating is still important at small stellar distances 
and near the midplane. The temperature in such disks is computed under the assumption 
that the upper surface of the disk is perfectly concave and smooth at all radii, which is a 
very good description of an unperturbed disk because thermal and gravitational instabilities 
are damped very efficiently (D'Alessio et al. 1999b). 

We are interested in planet formation inside such disks, which means that we might be 
forced to abandon the above assumption of a perfectly smooth surface. The newly formed 
planet core will distort it, affect the heating and cooling of the disk locally, and could have 
significant consequences for the further growth and migration of protoplanets. As discussed 
in Sasselov & Lecar (2000), the distortion need only be large enough compared to the grazing 
angle at which the starlight strikes the disk. This small angle has a minimum at 0.4 AU 
and increases significantly only at very large distances. The depth of the depression due to 
the additional mass of the planet, m p , will be proportional to (rn/h) 3 , where rn is the Hill 
radius, and h is the local scale height, with a shade area dependent on the grazing angle 
(Sasselov & Lecar 2000). 

Generally planet formation has been treated numerically in 2-dimensional disks; to 
date, there is only one recent study of the 3-dimensional effect due to planets embedded in 
a protoplanetary disk (Bate et al. 2003). Bate et al. model the 3-dimensional response of a 
gaseous viscous disk with a simple locally isothermal equation of state, without taking into 
account irradiation from the central star. 

In this paper we study in detail the development of a compression in a standard passive 
disk and the radiative transfer in it. We consider this to be the first step in developing a 
global 3-dimensional simulation, similar to the one by Bate et al. (2003), but with a realistic 
equation of state, especially near and inside the protoplanet's Roche lobe. We are explicitly 
interested in the planetary growth process before a gap in the disk is formed, hence our 
approach is limited to m p < 30 M ffi , which results in a small perturbation to the irradiated 
surfaces of a passive disk. 
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2. Heating Sources 

A schematic of our adopted disk model is shown in Figure 1. The disk is flared, so that 
disk height increases faster than distance to the star, consistent with observational evidence 
and theoretical models (Kenyon & Hartmann 1987; Chiang & Goldreich 1997; D'Alessio 
et al. 1998). We assume that the dust is well mixed with the gas and is the primary source 
of opacity at wavelengths characteristic of stellar emission. The dust absorbs incident stellar 
radiation at the surface and re-emits at longer wavelength, where the opacity is much lower. 
We are interested in the region of the disk that is optically thick, at distances of < 10 AU 
(D'Alessio et al. 1999a). This gives rise to three different layers with smooth transitions 
in the disk: the uppermost is optically thin to both stellar and disk radiation, the interior 
is optically thick to both stellar and disk radiation, and the middle layer is a transition 
region that is optically thick to stellar radiation, but still optically thin to disk radiation. 
We shall refer to this optically thick/thin region as the photosphere. The reprocessing of 
stellar radiation in this region is important for determining the interior temperature. 

Chiang & Goldreich (1997) calculate the temperature of only the thin/thin and thick/thick 
regions using energy conservation considerations, without detailing the temperature struc- 
ture of the photosphere. D'Alessio et al. (1998,1999a) calculate the detailed temperature 
structure of the photosphere in disks around young stars in a self-consistent way, includ- 
ing radiative transport and convection. They find that radiative transfer is the dominant 
mechanism for energy transport in the disk photosphere, and the primary heating sources 
are stellar irradiation and viscous heating, as shown by comparison to Calvet et al. (1991), 
who use the Milne-Strittematter treatment of the superposition of solutions (Milne 1930; 
Strittmatter 1974). We shall consider stellar radiation and viscous heating to be the pri- 
mary sources of heating, where the total temperature is arrived at by calculating the sum of 
the viscous and radiative fluxes, as T 4 = T* + T r 4 , where T v is the viscous temperature and 
T r is stellar radiation heating temperature. 



2.1. Viscous Heating 

The temperature due to viscous heating for a constant flux, gray atmosphere is 

T, 4 = !^(r, + 2/3) (1) 

where Ob is the Stefan-Boltzmann constant, is the optical depth of the disk's own radiation, 
and the viscous flux F v at a distance a, for a star of effective temperature T*, mass M*, and 
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radius i?*, accreting at a rate M a is 



3GM*M ( 
47ra 3 



( 



i2* 



) 



1/2' 



F v = 



a 



1 



(2) 



a 



(Pringle 1981). 



2.2. 



Stellar Irradiation 



We shall use a modified form of the Milne-Strittmatter treatment to calculate the effect 
of heating by stellar irradiation on a three-dimensional perturbation on a passive accretion 
disk. We assume that the dust is well-mixed with the gas in the disk and that they are 
thermodynamically coupled so that the dust temperature and gas temperature are the same. 
In general, there is no guarantee that these two temperatures are the same, or that the dust 
is either well mixed or uniform in composition. 

Following the Milne-Strittmatter treatment, the incident radiation from the star is con- 
sidered to be at some characteristic wavelength, typically in the visible, and the radiation 
emitted by the disk will be at a longer characteristic wavelength, typically in the infrared. 
The incident, short wavelength radiation will be indicated by the subscript s (stellar or scat- 
tered), and the re-emitted, long wavelength radiation will be indicated by the subscript d 
(disk or diffuse). Some fraction of the energy of the incident stellar radiation, a, is scattered 
at the same frequency. The remaining fraction, a — 1 — a, is absorbed by dust in the disk 
and re-emitted at a longer characteristic wavelength, typically in the infrared. This gives 
the equation of radiative transfer for the stellar short wavelength radiation, 



where p is the density, k s is the opacity at short wavelengths, and E is the incident energy 
flux given by 



— k • VI 8 = I a - 



aE e Ts a J I s dVt 



(3) 



K s p 





(4) 



The zeroth moment of this equation is 

—V • F s = aAnJs - (xE e- T % 
n s p 

and the first moment, making use of the Eddington approximation, is 



(5) 




(6) 
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The radiative transfer equation for the radiation absorbed and re-emitted within the 
disk is 

k-VI d = I d -B (7) 

where K d is the optical depth at disk radiation wavelengths. The zeroth and first moments 
(using the Eddington approximation) are 

— V-F d = Air(J d -B), (8) 

K d p 

— VJ d = 3 -^. (9) 

n d p 4tt 

2.2.1. Plane-parallel Disk 

In a plane parallel disk atmosphere, this problem is reduced to one dimension, and we 

can express everything as a function of optical depth normal to the surface, r d . We can 

substitute ^ 

V -> — and F M -> F s4 . 

n d p dr d 

If the angle of incidence of the stellar radiation is cos -1 p , then 

r d = T s p /q (10) 

where q is the ratio of opacities 

The condition of zero net flux relates the short and long wavelength fluxes as 

F s + F d = E p exp(-qr d /p ) (12) 
and the boundary conditions at the surface are that of an isotropic radiation field, 

2ttJ s (0) =F s (0) and 2vr J d (0) = F d (0). (13) 

Now, equations (5), (6), (8), (9), and (12) are a closed system of linear differential 
equations that we can solve directly for the temperature, T r = {i\B j ob) 1 ^ ■ The solution is 
given in Calvet et al. (1991) as 

B = ^^[C[ + C> 2 exp(-qr d /p ) + C 3 exp(-/3gr,)] (14) 



- 6- 



where (3 = y3a, 



and 



C[ = (l + C 1 )[2 + ^)+C 2 (2 + ^), (15) 



q J V Pv 



,2 



^ = ^(g-^), (17) 



a{2 + 3/i ) 

0(1 + 2/3/3) (1-/3 W 1 ' 



Figure 2 shows how temperature varies with fj, at the surface (r d = 0) and in the interior 
(t<2 — > oo). The temperatures are normalized to an effective radiant temperature T = 
(i^o/cs) 1 ^ 4 ! which is the blackbody temperature corresponding to the incident energy flux. 
Note that at the surface the temperature remains relatively constant with /io and is greater 
than T due to the absorptive properties of the dust, in accordance with the "superheated" 
surface layer proposed in Chiang & Goldreich (1997). In contrast, the interior temperature 
is sensitive to the incident angle, especially at grazing incidence. Thus, perturbations in the 
surface that change the angle of incidence can significantly affect the interior temperature. 

By definition, = 2/3 is where the disk becomes optically thick to its own radiation. For 
q > 1, the temperature at = 2/3 is can be approximated by the interior temperature, since 
as qr d becomes large, the exponentials in equation (14) vanish and B becomes independent 
of Td- This indicates that in the absence of viscous heating, the disk interior is vertically 
isothermal for r d > 2/3. 



2.2.2. Disk with Perturbation 



Now we consider radiative transfer on a perturbed disk. A point P(x, y, z) within 
the disk "sees" radiative flux coming from the surface of the disk. Each area element 5 A 
contributes a solid angle SO, of flux to P. Figure 3 shows a schematic of such an area element. 
The angle of incidence of radiation on this surface element is cos^ 1 /i. If k is the vector from 
5 A to P and n is the unit normal to 5 A, then v is the cos of the angle between n and k. 
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We approximate the flux contribution from 8 A as if it were a surface emitting isotropi- 
cally with intensity 

= F d (T d ,fl) 
7T 

and the corresponding contribution to the total flux is d¥ = Ik dfl. Then the total flux at 
P over all surface contributions is 

F tot = - [ F d k5Q. (20) 



7T 



Since V is a linear operator, 

1 fdF d 



-V • F tot = - 



Defining 



h-k5tt = 4n(J- J J d v5Q-^- J Bv5t?J . (21) 



J tot = i y J d (r d , n) v SQ and (22) 

But = B(T d ,n)v5Q (23) 

we see that F to t, <Aot, and -B to t satisfy equations (8) and (9), and we can calculate the 
perturbed temperature as 

Tr= ^^y /4 . (24) 

For a plane-parallel surface, r d and /i are constant, and v = cos 9, where 9 is the angle 
with respect to the surface normal. We then recover the solution for the plane-parallel case. 

This method of calculating radiative transfer can be applied to any three-dimensional 
disk configuration, including perturbations induced by protoplanets. Sufficiently massive 
protoplanets will open a gap in the disk, in which case the effects of shadowing and illumi- 
nation across the gap can be calculated for all azimuthal angles. In the following section, we 
will calculate radiative transfer on small perturbations induced by planets too small to open 
a gap, but large enough to act on the disk locally. 



3. Hydrostatic Equilibrium 

Now we consider a perturbation induced by the gravitational potential of a protoplanet 
on a circumstellar disk. We shall consider only the gravitational effects of the protoplanet 
on hydrostatic equilibrium and assume that any resonant effects are washed out by gas drag. 



- 8- 



3.1. Vertical Density Profile 

In order to find the vertical structure of a non-self-gravitating gaseous disk orbiting a 
central star, we solve the equation of hydrostatic equilibrium: 

VP = -pV$ (25) 

where P is the gas pressure of the disk and $ is the gravitational potential of the central 
star. Since we are interested in the vertical structure, we consider only the z components 
and assume that the disk is isothermal in the z direction. For a disk without a perturbing 
protoplanet, $ = GM+/a where M* and a are the mass and distance of the star, respectively. 
If c s is the isothermal sound speed, then 

2 d\np GM* 

C ° — = -— Z - (26) 
If we assume that a^> z and define h = c s a 3 / 2 / (GM*) 1 / 2 = (c s /t>^)a, then 

p(z) = p exp (-^2) > ( 27 ) 

i.e. the density has a Gaussian distribution with a scale height of h which is determined by 
the local temperature. 

The insertion of a planet into a passive disk adds an additional term to $, representing 
the potential due to the planet: 

$ = GM^/a + Gm p /r p 

where m p and r p are the mass and distance of the planet. Assuming that the planet is in 
the midplane (z — 0), Eq. (26) becomes 

d mp z ( M p a 3 \ z ^ 



dz h 2 \ M+h 2 J rl 



The Hill radius is defined as r H = (M p /3M^) 1 / 3 a, and this equation integrates to 



p(z) = pi exp 



z 2 3r 3 
+ ' 



2h 2 h 2 r p 



(29) 



where r p = ^ x 2 + y 2 + z 2 , with the planet as the coordinate origin. So long as r H h and 
\z\ > 7"h, the effect of the planet is a small perturbation on the density profile. 
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We normalize this equation by matching the density at z = to the unperturbed density, 



thus 



P P ( Z ) = Poexp 



z 2 3r 2 
+ " 



2h 2 h 2 V + y 2 + z 2 y 1 x 2 + y 2 



(30) 



Note that there is a singularity ar r p = 0. Figure 4 shows the shape of this density profile. 
In Figure 4a, we hold r H //i fixed at 0.5 and vary y 'x 2 + y 2 . For a disk with ft, = 0.04 AU at 
1 AU around a central star of 0.5 M , this corresponds to a protoplanetary mass of 4 M e . 
For [x 2 + y 2 ) > r^, the perturbed density profile does not deviate significantly from the 
unperturbed density profile. For (x 2 + y 2 ) < r^, the density profile develops a sharper profile 
at z < h. In Figure 4b, we hold a/x 2 + y 2 fixed at the Hill radius and vary the mass of 
the protoplanet, effectively varying rn- Larger protoplanets have the effect of decreasing the 
overall density profile, whereas position controls the shape of the profile, particularly close to 
the protoplanet. Since the flow of gas within the Roche lobe cannot be adequately described 
hydrostatically, we exclude the region ^ x 2 + y 2 < th from our analysis. 

The normalization on equation (30) is somewhat arbitrary, but at the disk heights that 
are of interest to us, the overall normalization has small effect. We demonstrate this as 
follows: as z > r H , equation (30) becomes 



P P {z > r H ) ~ p exp 



34 



2h 2 h 2 ^x 2 + y 2 



(31) 



that is, a Gaussian with respect to z with some normalization that depends on r H and 
^/x 2 + y 2 . The height of the photosphere, H, is generally several times the pressure scale 
height. If we take H/h = 5, then the change in density with vertical distance is 



Ap 
P 



Az dp 
p dz 



(z = H) 



25Az 
H ' 



(32) 



So, a density change by a factor of 1/2 corresponds to a shift in vertical height of 2%. In other 
words, the density gradient in the photosphere is so high that changes in the normalization 
of the density will not significantly affect our results. 



3.2. Photosphere Height 

We calculate the height of the perturbed photosphere from the perturbed density profile 
as the isodensity contour at the density of the unperturbed photosphere. That is, if H is 
the unperturbed photosphere height and H = H (l — e) is the new photosphere height, then 
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we set 



or 



Pd (H ) = Pp {H) 

m h 2 , ( i 

2 2 H I + y i + H 2 



\/x 2 + y 2 



(33) 



Setting ^ x 2 + y 2 = 7-f/o, then to first order in e 



H « # o < 1 



(34) 



so that the depth of the perturbation scales as (rn/i^o) 3 - For typical disk models, the ratio 
of the photosphere height to the thermal scale height, H /h, is between 3 and 5, depending 
on the distance from the star and is often taken to be a constant (Kenyon & Hartmann 1987; 
Chiang & Goldreich 1997; D'Alessio et al. 1998). Restating the scaling of the depth of the 
perturbation as oc (ru/h) 3 , we see that this agrees with the estimate made in Sasselov & 
Lecar (2000). 

The shape of this perturbed photosphere is shown in Figure 5 for varying values of 
r H . The perturbation is small for r > r H , and there is a singularity at r = 0. As the Hill 
radius increases, the depth of the perturbation also increases. Note that if H /h — 5, then 
r-a/Ho = 0.2 means that the Hill radius is comparable to the disk pressure scale height and 
the perturbation is no longer small. Protoplanets this massive are likely to open a gap in 
the disk rather than simply inducing a perturbation to the photosphere. 



4. Disk Shear 

The differential rotation of the disk causes the disk material near the protoplanet to 
move through the perturbation at a rate that depends on the velocity with respect to the 
protoplanet's orbit. Material moving along a given streamline in the disk may pass through 
the perturbation too quickly to experience significant heating and cooling, thereby dimin- 
ishing the effect of the perturbation on the temperature structure of the disk. 

The velocity of the disk material with respect to the protoplanet is v ~ a(Q — Q p ) 
where Q and Q p are the orbital angular velocities of the disk material and the protoplanet, 
respectively, and a is the orbital radius of the protoplanet. The orbital angular velocity of a 
gaseous disk is 

r p dr 



- 11 - 



In thin disks, such as we are investigating, the pressure gradient term is small so the orbital 
velocity is nearly equal to the Keplerian velocity. 

The heating/cooling rate of the disk material can be expressed as 

FfT 

C^- = F-aT 4 (36) 
dt 

where we take F = u{T 4 + T*) and C is the specific heat per unit surface area of the disk. 
We shall adopt a specific heat of C = kH/m where k is the Boltzmann constant, S is the 
total surface density, and m is the mean molecular mass. 

We shall assume that the streamlines of the disk material follow equation (35), and 
that they are not perturbed by the protoplanet's graviational potential. This is an adequate 
approximation for disk material outside the protoplanet's Hill radius. Along each streamline, 
we calculate the total radiative flux, F, at a given position, and using the velocity of the 
streamline with respect to the protoplanet along with heating/cooling rate, we can calculate 
the steady state temperature at each position in the disk. 



5. Model Parameters 



Dust properties can be parametrized by a, k s , and q. Then, the temperature profile in 
a plane-parallel disk is determined completely by F v , E , /x , and r d , where r d is measured 
normal to the surface of the photosphere. Stellar properties determine F v and E , and the 
disk profile determines /xo- If Hq oc r^, then the incident angle is determined by 

a = ^ ~ l)Ho/a (Z7) 

^ [l + (H /ani + e(H /a)i] 1 1 



A perturbation changes the local angle of incidence and optical depth. In particular, 
without plane parallel symmetry, there is no single optical depth that parametrizes the 
distance below the surface. Instead, we sum over different lines of sight, using the optical 
depth along each line of sight to determine the contribution to the disk flux as in equation 
(23). The local angle of incidence depends only on the geometry of the perturbed surface, 
while the optical depth also depends on the density structure along the line of sight. The 
shape of the perturbed surface depends on r H and H , and the density structure depends on 
the disk pressure scale height, h. 



- 12 - 



5.1. Fiducial Model 

Motivated by observations of T Tauri stars, we will assume for our fiducial model that the 
central star has mass M* = 0.5 M Q , radius i?* = 2R Q , and effective temperature T* = 4000 
K, and that the protoplanet is at a = 1 AU from the star. We assume an accretion rate of 
M a = 10~ 8 M yr -1 , so that T v = 72 K at r d = 2/3 at 1 AU. With these parameters, heating 
in the photosphere will be dominated by stellar irradiation rather than viscous heating. 
Near the midplane, however, viscous heating will dominate. Since we only consider the 
photosphere in our model, the calculation of the midplane temperature is outside the scope 
of this paper. 

We will assume that the fraction of scattered radiation is a = 0.28, with opacity in the 
optical of k s = 400 cm 2 g _1 , and ratio of opacities q = 20. 

We shall assume that H = 0.16 AU(r/l AU) 9/7 , consistent with e.g. D'Alessio et al. 
(1999a). In calculating the detailed density structure in the photosphere, we shall assume 
that h — 3.4 x 10~ 2 AU so that H /h — 5. A protoplanet with mass m p — 1 M ffi (10 M ffi ) 
will have a Hill radius of 1.26 x 10~ 2 AU (2.72 x 10~ 2 AU). As shown by Bate et al. (2003), 
planets with masses < 30 M e are not massive enough to open a gap, and our perturbative 
treatment is valid. We have calculated models with lM e , 10 M e , and 20 M ffi planets, and 
present the results below. 

5.2. Numerical details 

We define our coordinate axes so that x is aligned with the radial direction, y is aligned 
with (i.e. the direction of planet's orbit), and z is perpendicular to the orbital plane. The 
coordinate origin is set at the surface of the disk above the planet position, so that the 
planet's coordinates are (0,0, —H ). 

We calculate the shape of the disk surface in the x and y directions over a range — 10rn < 
x < 10r H and — 10r H < y < 10r H , and use this surface to numerically integrate temperatures 
within the photosphere, to a depth of zir^ = 2/3) = —0.011 AU. The temperature is 
calculated at each point within a 64 x 64 x 64 grid. 

When summing over surface elements to find B tot as in equation (23), we can approx- 
imate 5Q ~ u5A/£ 2 for 5 A <ti £ 2 , where £ is the distance between the point P and 5 A. 
However, for points close to the surface, this approximation breaks down. In addition, v 
may have large excursions over 5 A. For this reason, we assume that each 5 A is approxi- 
mately planar, so that B{jd, /i) does not vary, but we calculate f sn v d£l analytically from 
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the limits of 5 A. In other words, we calculate 

B tot = B(T d ,fi) Jvdn (38) 

and solve for the temperature. 



sn 



6. Results 

Figures 6, 7, and 8 show the spatial distribution of temperatures in the photosphere of 
the fiducial model for 1 M e 10 M e , and 20 M e planets, respectively. The colors represent 
temperatures in Kelvin, as indicated by the colorbars. The horizontal axis indicates increas- 
ing radial distance from the star so that the photosphere is illuminated from the left. The 
cross-section is taken at the Td = 2/3 surface. The bottommost row and leftmost column 
show the unperturbed temperature of 133.6 K for comparison. This temperature is also 
indicated in the colorbar by a dotted white line. 

The direction of the planet's orbit is upward on the plots, so the interior disk material 
moves faster, and the exterior disk material moves slower than the planet. This causes the 
heated and cooled regions to shear out. The positions of minimum and maximum temper- 
atures are offset from the planet's position since it takes some time for the heated/cooled 
material to come back to the equilibrium temperature. 

For a 1 M e protoplanet, the shadowing and illumination at the surface results in a 
slight temperature change, with the area to the right of the planet reaching a maximum 
temperature of 135.5 K, and the area to the left cooled to a minimun of 130.7 K - a change 
of —2.2% to +1.4%. Thus, we see that an Earth-mass planet is too small significantly affect 
the interior temperature of the disk. 

A larger planet will create a larger perturbation, which will lead to a greater temperature 
change. Figure 7 and Figure 8 show the temperature distribution in the photosphere for a 
10 M ffi planet and 20 M ffi planet, respectively. The plots are scaled to the Hill radius, which 

1 /3 

goes as m p , so the physical range in Figures 7 and 8 is actually larger than that in Figure 
6. The high temperatures close to the Hill radius are due to the low density in this region, 
which gives a low optical depth to the illuminated surface. However, since this region is 
not in hydrostatic equilibrium, we exclude it from our analysis. Outside the Hill radius, the 
temperatures at r d = 2/3 for a 10 M ffi protoplanet range from 123.7 K to 142.3 K, or by 
—7.4% to +6.5%. For a 20 M e protoplanet, which is close to but not at the limit where a 
gap in the disk can form, the temperature range increase slightly, with a minimum of 122.0 K 
and maximum of 146.6 K, or by —9.7% to +8.7%. Thus, a protoplanet with mass 10 M ffi or 
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larger may have a significant effect on the local temperature structure of a protoplanetary 
disk. 



7. Discussion 

As planet-forming disks are likely to be passive (Calvet et al. 2002), we have explored 
stellar irradiation effects due to a disk surface perturbation induced by an embedded proto- 
planet. We find that even a relatively small protoplanet (in the range from 1 to 20 Earth 
mass) produces enough of a perturbation to cause shadowing. In particular, we find that in a 
standard passive disk accreting at lO~ 8 M yr _1 , the temperature variation induced by a lM e 
planet at 1 AU is 5 K, or only about 4%. However, we find that even though a protoplanet 
may not be massive enough to open a gap in a disk, it may still be able to significantly 
change the temperature structure in a the disk material near it. A 10 M e (20 M e ) planet 
can induce a temperature variation of 19 K (25 K), or 14% (18%). 

This results in a local temperature perturbation which is large enough to affect the 
local properties of the disk and the accretion rate, and possibly the migration rate of the 
protoplanet. Although the perturbation in temperature may be too subtle to directly observe, 
it can have serious consequences for planet growth and migration. 

Ward (1997) has shown that the rate of Type I migration is strongly dependent upon 
the temperature gradient of the disk. If the temperature decreases with distance from the 
star, then the total net torques result in inward migration of the secondary. However, if 
k = —dlnT/dlnr < —1, then the opposite is true and the secondary will migrate outward. 
The temperature perturbation near a protoplanet locally decreases the value of k, and may 
result in halting or even reversing Type I migration. This may resolve the conundrum that 
Type I migration timescales are typically much less than observed disk lifetimes of ~ 10 7 
years (Hartmann et al. 1998). Yet to date, no simulation of planet-disk interactions accounts 
for local temperature pertubations as described here. 

We have not yet considered the behavior of the dust (as a function of size) in the 
perturbed region. We expect the dust to remain coupled to the gas and compress its scale 
height accordingly. In fact crude estimates of dust sedimentation in the vicinity of the planet 
confirm that. However, a detailed estimate of the sedimentation timescale is necessary in 
order to evaluate the amount of residual dust, e.g. swept by stellar radiation from the inner 
edge, and its total opacity. These are issues we plan to address in future work. 

In a forthcoming paper, we will undertake a parameter study of our model for radiative 
transfer on a perturbed disk. By varying attributes of the disk, star, and protoplanet, we 
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will investigate the possibility that a relatively small protoplanet which is not large enough 
to open a gap in a disk might still be able to affect disk structure. 

We thank Nuria Calvet for numerous helpful comments in the preparation of this paper. 
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Fig. 1. — Schematic of a flared disk. Proportions are not to scale. The labels thick and thin 
refer to optical depth in short and long wavelengths. 



-17- 




Fig. 2. — The temperature due solely to stellar irradiation at the surface of the disk (dashed 
line) and in the interior (solid line) for varying values of the angle of incidence, /j,q, where 
T = (£ /^) 1/4 . 
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Fig. 4. — Perturbed vertical density profile for varying parameters. Density is normalized to 
the midplane density. Both plots are to the same scale and show the unperturbed density 
profile as a solid line, (a) Density profile for fixed Hill radius and varying values of (x 2 +y 2 ) 1 ^ 2 
as indicated, (b) Density profile for (x 2 +y 2 ) 1 ^ 2 equal to the Hill radius and varying planetary 
masses. 
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Fig. 5. — Perturbed photosphere height as a function of distance in the xy-plane, r, for 
varying Hill radii. The solid vertical line indicates where r = th- 
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Fig. 6. — Temperature profile of the disk photosphere near a protoplanet of mass 1 M ffi at 
Td = 2/3 in the xy-plane. The colors indicate the temperature in Kelvin, scaled as shown in 
the colorbar. The cross in the center represents planet position, with the arrow indicating 
the direction of the planet's orbit. The blacked-out circle indicates the extent of the Hill 
radius. The dotted white line in the color bar indicates the unperturbed temperature at 
r d = 2/3. 
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Fig. 7. — Same as Figure 6, for a protoplanet of mass 10 M ffi . Note that the dimensions on 
the figures scale with the Hill radius, so the physical sizes of the plots are 10 1 / 3 larger than 
in Figure 6. 
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